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Abstract. General expression for the thermodynamic potential of the model of semi¬ 
infinite jellium is obtained. By using this expression, the surface energy for infinite 
barrier model is calculated. The behavior of the surface energy and of chemical 
potential as functions of the Wigner-Seitz radius and the influence of the Coulomb 
interaction between electrons on the calculated values is studied. It is shown that 
taking into account the Coulomb interaction between electrons leads to growth of 
the surface energy. The surface energy is positivein the entire area of the Wigner-Seitz 
radius. It is shown that taking into account the Coulomb interaction between electrons 
leads to a decrease of the chemical potential. 
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1. Introduction 

The development of quantum-statistical theory of Fermi systems with interfaces is one 
of the most important problems of contemporary statistical physics. In particular, the 
richness of surface phenomena and the rapid development of experimental methods of 
investigation of surfaces requires the development of theory of such systems. 

The most popular theoretical method for studying in the area of research is the 
density functional theory m. that have been developed from the well-known Thomas- 
Fermi method for atoms. By construction, the density functional theory is the one- 
particle approach and can not properly take into account the many-body correlation 
effects. Therefore, the energy functionals for inhomogeneous systems are mostly used 
in the local density approximation [2], namely the electron density distribution n(r) is 
substituted by the average electron density n = const in the well-known expressions 
of the theory for homogeneous systems. This approach is questionable |1], since the 
presence of the interface brings both quantitative and qualitative changes of various 
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characteristics of an electronic system, e.g. the image forces, that cannot be obtained 
from the density fnnctional theory in principle. 

Density fnnctional theory has a characteristic problem of snrface energy. Namely, 
the snrface energy of semi-infinite jellinm calcnlated in this theory tnrns ont to be 
negative for large valnes of the electron concentration (rg < 2.5 ob, where Vg is the 
Wigner-Seitz radins) |5]. This is physically incorrect. The snrface energy mnst be 
positive, otherwise the metal wonld spontaneonsly split. At present, general belief is that 
the canse of negativity of the snrface energy is the replacement of the discrete ionic lattice 
by nniform positive backgronnd. Thus, in the work of Lang and Kohn [5] a discrete 
lattice is accounted for using the first-order perturbation theory in the pseudopotential. 
As a result, the surface energy becomes positive and is satisfactorily consistent with the 
experimental data for a number of simple metals j5]. A variational procedure has been 
developed by Monnier and Per dew PIT] to take into account the averaged effect with 
the introduction of discrete additions to the potential inside a metal, which depends 
on the structure of a lattice and the surface. The obtained results for the surface 
energy were very close to the results of Lang and Kohn and were in better agreement 
with the experimental data. Later, Appelbaum and Hamann [8] used the local density 
approximation and performed the calculations for the Cu(lll) surface, with the account 
of the discrete ionic lattice without perturbation theory. They have obtained good 
agreement with the experimental data. The authors of all these studies have assumed 
that the non-local exchange-correlation effects are negligible and can be omitted. In 
the works pen], the calculations of the surface energy were performed with varying 
position of the last (exposed) ion layer and good agreement with experimental data 
for the surface energy of simple metals was obtained. In the works [IlHl3], a stabilized 
jellium model has been proposed. In this model the pseudopotential correction, which is 
averaged over the Wigner-Seitz cell, is incorporated into the effective potential inside the 
metal. This model yields positive values for the surface energy. Thus, the consideration 
of discrete ionic lattice permits to solve the problem of negative values of the surface 
energy. 

However, it is not clear that neglecting the discreteness is the principal and only 
reason of discrepancy between the theory and experiment. Probably the theory could 
be improved, still remaining in the framework of the jellium model. There have been 
attempts to go beyond the local density approximation. Namely, Schmit and Lucas [13], 
Craig na and Peuckert [I6] have considered certain non-local contributions to the 
exchange-correlation term to the surface energy, due to the change of zero energy 
of plasmons and the appearance of surface modes during separation of crystal into 
fragments. From a good agreement of this contribution with the experimental data for 
the surface energy, the authors of these works proposed to identify this contribution 
with the total surface energy, suggesting that other contributions, that have not been 
accounted for, canceled each other. This approach is actively debated in Refs. UMl]. 
The expressions for the exchange-correlation energy of bounded electron gas have been 
obtained by Harris and Jones [22l[23], Wikborg and Inglesfield [25], Johnson and 
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Srinivasan [27], focusing in the analysis of non-local effects. The exchange part of 
the surface energy was calculated for electrons in a potential box with inhnitely high 
walls [221123], and the exchange-correlation and exchange parts were calculated in the 
random phase approximation [25] . A comparison of these results with those calculated in 
the local density approximation has showed that the exchange-correlation parts differ by 
about 10%, the exchange parts differ by 50%, and the correlation parts differ by 6 times. 
The local density approximation works much better for the sum of the exchange and 
correlation parts of energy rather than for each individual contribution (see also IZSD- 
Further calculations, using more realistic models of the surface barrier [281 - I32] . namely 
the gradient expansions [321139] . the analysis of Langreth and Perdew [381139] . have found 
a decisive contribution of the local density approximation into the exchange-correlation 
part of energy as well. In particular, the calculations [35] show that the non-local 
corrections do not exceed 16%, although the relative contribution of non-locality into 
the total surface energy can be much larger (up to 40%) because for many metals 
the exchange-correlation part of the surface energy is greater than the total surface 
energy [5]. 

More recent studies have renewed the debate about the correctness of the 
application of the local density approximation in the calculation of the surface energy. 
Using the Fermi hypernetted-chain equations [iQlST] have obtained signihcantly higher 
values of the surface energy than calculated by Lang and Kohn in the local density 
approximation [5]. In contrast, the values of the surface energy obtained in the 
calculations using density functional theory with nonlocal functional of Langreth 
and Mehl [43] , are much closer to the results obtained in the local density approximation. 
Calculations of the surface energy by using quantum Monte Carlo method [441145] 
have shown that the values of the surface energy obtained at high concentrations 
{ts ^ 2.07 Ob) are in a good agreement with the results obtained by using nonlocal 
functional, but at lower concentrations Cg ^ 3.25 ob) they are in a good agreement with 
the results of Krotscheck and Kohn [401I4T] . On the other hand, Pitarke has concluded 
that the local density approximation leads to a small error in the exchange-correlation of 
the surface energy [46]. In this work, the long-range correlations are taken into account 
self-consistently in the random phase approximation, while the short-range correlations 
are included in the time-dependent local density approximation [471148] . 

In Ref. [49] the surface energy is calculated by using the one- and two-particle 
distribution functions of electrons, that are obtained in [50] taking into account the 
Coulomb interaction between electrons. It is found that the surface energy is positive 
in the entire region of electron concentration. 

Takahashi and Onzawa have calculated the electron density distribution and the 
surface energy of non-interacting electron gas for hnite barrier model [51]. It is 
interesting that the electron density without self-consistency is very close to the self- 
consistent calculations of Lang and Kohn [5], and the surface energy is positive for all 
concentrations of electrons. Moreover their result is close to the result of Lang and 
Kohn at low concentrations. 
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In the present work, our principle objective is to construct a consistent quantum- 
statistical theory of a simple metal with the interface ‘metal-vacuum’ in the framework 
of the jellium model. An effective potential of inter-electron interaction, which we 
have studied recently in [521156] . is crucial for calculating of the general expression 
for thermodynamic potential. It is shown that within certain approximations the 
thermodynamic potential can be represented as a functional of the one- and two- 
particle distribution functions of electrons. At low temperatures, the nonlinear algebraic 
equation for the chemical potential and a general expression for the internal energy 
are obtained from the thermodynamic potential. It is shown that obtained equation 
for the chemical potential of non-homogeneous system is similar to the equation for a 
homogeneous system. The chemical potential is calculated as a function of the Wigner- 
Seitz radius. The one- and two-particle distribution function of electrons are calculated 
for the inhnite barrier model. The expressions for the extensive and surface contributions 
to the internal energy are obtained. For the same model, the calculation of surface 
contribution to the internal energy, which is the surface energy at low temperature, 
is performed. The behaviour of the surface energy as a function of the Wigner-Seitz 
radius is studied. We have found that the surface energy calculated in the present work 
is positive in the entire concentration range typical for metals, and at low concentrations 
is consistent with the calculations of Lang and Kohn [5]. The influence of the Coulomb 
interaction between electrons on the calculated characteristics are studied. Detailed 
calculations of these properties is given in Ref. IS]. 


2. Model 


We consider a system of N electrons in the volume V = SL in the held of positive 
charge with the distribution 

ad,(R||,Z) = = eoB(-Z) = I ^ J ° , (2,1) 

where Z = 0 is the dividing plane, 9{x) is the Heaviside step function, R|| = {X,Y), 
X,Y G [—y/S/2,+y/S/2], Z G [—L/2,+L/2]. The condition of electroneutrality is 
satished. 


-\-Lj2 


lim 

S^L —^cx) 


dR|| /dZ ^jeu(R||, Z) = eN, e> 0, 


S -L/2 

moreover in the thermodynamic limit we have, 

eN 


lim —— = lim 

N,S,L^oo SL N,V^oo V/2 


— Qo- 


( 2 . 2 ) 


(2.3) 


This model system is known as “semi-inhnite jellium” and it is one of the 
simplest models of semi-inhnite metal, which satisfactorily describes simple metals. The 
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Hamiltonian of the model is, 


Hjell - 


2m 


N 


N 




2=1 


2 ^ |r,' 

i^j=i 


N 

^/dR 

j=^v 




Rl 


2j 7 IR 1 -R 2 I 

V V 


(2.4) 


where is the position of j-th electron; the hrst term is the kinetic energy of electrons 
(m is the electron mass), the second term is the potential energy of the inter-electron 
interaction, the third term is energy of interaction of electrons with the positive charge, 
the fourth term is potential energy of the positive charge. 

From the Hamiltonian fl2.4p we extract a Hamiltonian of the inhnite jellium 
model Rjeu^ 


N 


N 


rrunif 
^jell - 


2m, 2 Irj — r,j 

i^j=i i=i 


dR 


2=1 


V 


e^N/V 
ki - H.| 


27 J "“IRi-RjI’ 

V V 

here physical meaning of the terms are similar to the terms of the Hamiltonian 
Thus we get. 


(2.5) 


N 

V V 


^jell(Rl)^'jell(R2) — {eN/Vy 


|Ri ~ R 2 


where 


^surf(rj) = J 

V 


e{eN/V - Pjeu(R)) 


ki - R-l 


( 2 . 6 ) 


(2.7) 


is the surface potential acting on the electron. This potential is formed by the deviation 
of the positive charge distribution from the uniform one. So, if instead of pjeii(R) we 
put the uniform distribution eN/V, then the surface potential and the last term in the 
Hamiltonian (12.bh disappear and one obtains. 


1™.. ^jell = 


It should be noted that as a consequence of the symmetry of the model, the surface 
potential Kurf(r) is a function of the normal to the dividing plane coordinates of the 
electron only, the motion of the electron in a plane parallel to the dividing plane is free. 


i.e.. 


f4urf(r) = Kurf(2^)- 


In order to calculate the thermodynamic potential of the system, it is convenient 
to present the Hamiltonian (12.61) in the secondary quantization representation. 
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3. Secondary quantization representation 


We introduce the single-particle wave functions 'ha(r) and the corresponding energies 
Ea of the electron in the held of the surface potential 14 urf( 2 :), 


2m 


^ + i4urf(2^) 


^,(r) = 


(3,1) 


which we use to construct the representation of the secondary quantization. 

Since the potential in the stationary Schrodinger equation fl3.1|) depends only on the 
normal to the dividing plane coordinate of the electron, the variables can be separated. 
Then we obtain, 

hY 


Ea = 


2m 


+ a — (P; Ct), 






(3.2) 

(3.3) 


where r|| is two-dimensional coordinate of the electron in the plane parallel to the 
dividing plane, hp is the moment of the electron in this plane, and 

27rnr 


P = (Px,Py), Px,y = 


''X,y 


m 


= 0,±1,±2, 


(3.4) 


a is some quantum number, that depends on the form of the surface potential, the 
functions (pa{z) satisfy the one-dimensional stationary Schrodinger equation, 

d^ 


+ Kurf(^) 


(faiz) = ea(fa(z). 


2m dz'^ 

In the secondary quantization representation constructed by means of the wave 
, the Hamiltonian of the system becomes. 


functions 

H = ^«(p)4(p)«Jp) - Y, ^(q, + i"fc(q)pfc(q)p-fc(-q), (3.5) 


2SL 


p,a q^O q^O k 

where a]j(p), a„(p) are the operators of electron creation and annihilation, respectively, 
in the state (p, a), and the standard commutation relations are, 

{®ai (Pl) 1 (P 2 )} ,P2 *^ 01 ,02 ) (3.6) 




(3.7) 


p,a 


is the particle number operator, u (q, 0) = ^^^(q) = 47re^/(q^ -|- kf) is the Fourier- 

transform of the Coulomb interaction, = ^rUx^y, rOx^y = 0, ±1, ±2, ..k = 
n = 0, ±1, ± 2 ,..., 


pt(q) = (“I® ‘‘‘'l“')“!,(p)«„.(p - q) 


(3.8) 


p,a,a' 














Semi-infinite jellium: thermodynamic potential, chemical potential, surface energy 7 
is the mixed Fourier-representation of the local density of electrons, 

+ Z//2 

(3-9) 

-L/2 

It is worth noting that in the eqnation fl3.5p . there are no terms with q = 0, dne to 
the electronentrality condition fl2.2p . 

The Hamiltonian in the form (13.dp is convenient to calcnlate the thermodynamic 
potential by the fnnctional integration method. 


4. Thermodynamic potential 


4 . 1 . Functional representation 
The grand partition fnnction, 

S = Sp exp [ —/5(i4 —/iiV)], (4.1) 

that determines the thermodynamic potential of the system, 

n = (4.2) 

and other thermodynamic fnnctions, in the interaction representation becomes, 

S = Soexp (^^(Ar)o^i/(q,0))Sint, (4.3) 

q^O 

where Sq = Spexp ( — (3{Hq — /iiV)), Hq = i?Q,(p)a),,(p)aQ,(p) is the Hamiltonian of 

p,a 

non-interacting system, /i is the chemical potential, 

(...)o = 7^Sp(e-^(^°-^^)...), (4.4) 

{N)o = (4(p)a„(p))Q = Y na{p), 

p,Q: p,Q; 

na{p) _|_ I 

is the Fermi-Dirac distribntion, 

Sint = {S{fi))o, 


S{f3) = T exp 


2SL 


^ki(l)pk{q., /3')P-fc(-q, 

k 




(4.5) 

(4.6) 


T is the symbol of chronological ordering of “times” fi = 1/6, 6 is the thermodynamic 
temperatnre. 
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For further calculations it is convenient to switch to spectral representation, 


Pfc(q, = p j 


Pkiq.,/3') = 


where u = -^n {n = 0, ±1, ±2,...) are Bose frequencies. Then (14.51) becomes, 


(4.7) 

(4.8) 


S(P) = Texp 


1 


2SL 


EEE ''t(q)p/i(q. '')p-t(-q> -i') 


(4.9) 


qT^O k V 

In order to simplify the approximation iS(/9) according to fl4.4p we switch to 
functional representation for S{fi) [5HII59], using Stratonovich-Hubbard identity [60] . 

+ 00 


exp 


-jy'Ay 


= (detA)-^/2/ (dx)exp 


--x^A-^x + ix^y 
2 T 




(4.10) 


where, (dx) = ])([ y"'" = (|/i,... ,|/n), = (xi,... ,Xn), A is the positively dehned 

k=l 

matrix. Then, for fl4.9p we obtain. 


qT^O k V 

X /(dw)exp 




q^O k V 


X T exp ■EEE Wfc(q, i^)pfc(q, v) 

. qT^O k V 

where (den) is the element of the phase space, 

(du )=n n n 

q>0 k^O 


(4.11) 


TT 


TT 


ujk{ci, v) = a;fc(q, v) + in;fc(q, u), 

Wfc(q, = w^fc(-q,-z^), 

Wfc(q, 

Note, that due to the fact that operator variables pfc(q, v) are under sign T-ordering 
in fl4.1ip . it is impossible to perform integration by fi' in Eq. fl4.7p . 

By making average of S{(3) according to (14.4p . we obtain. 


■"int — 


{sw)o = 

q^O k V 


(4.12) 
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where, 




X exp 


E 

ri.^2 


J(a;)=exp - 2 

qT^O k V 

X] 9^fcl,...,fc„(qi,«^l,---,qn,Z^n)Wfel(qi,I^l)...X;fc„(qn,Z^n) , (4.13) 

Qn 7 ^ 0 , /Ctt, , Vji 

■ ■ ■ , qn, = i”(Tpfc,(qi, I/i) . . . Pfc„(qn, Z^n))o,c 

~ ^qi+q2 + ... + qn ,0 ^i/l+i/2 + ... + !2n,0 

are the so-called irreducible mean values (cumulants), 5 is the Kronecker delta, 
911^(q, z/) = 0, because q 7 0. 

In general terms, the calculation of integral fl4.12p is a complicated problem due 
to the exponential index having terms with u ^ 3. Their neglectance results in the 
Gaussian approximation (or the so-called random phase approximation). As a rule, 
calculation of this integral is done by means of series expansion of the non-Gaussian part 
of the integral fl4.13p with subsequent averaging with Gaussian distribution and partial 
summing up of the terms which give the most important contribution. In contrast, 
in [52] it was shown that the integrand Jipj) can be approximated by a Gaussian form 
Jq,{uS), introducing the unknown function i/)\ 

“ ^ z/))"V(q, i^)wfc2(-q, -v) 

q^O ki,k 2 i' 

where gkT_,k2i^-’ Fourier-transform of effective inter-electron interaction, 

(|^fci,fc2(q, = (Ai^fci(q))~^<^fci+fc2,o - ^fci,fc2(q, ^)- (4-i5) 

We will seek the unknown function Dki,k2{^.j from the condition that the mean 
value of a;fc^(q, z/)a;fc 2 (—q, —^^), calculated with the distribution Jiu), is equal to the 
mean value, obtained with Gaussian distribution Jq{u:), namely from the condition, 

a;fci(q, z^)cafc 2 (-q,-z/) = (ozfci(q, z^)cafc 2 (-q,-z/))^, (4.16) 

where we have introduced the following notation, 

f{duj)J{u)... 


Jg{uj) = exp 


(4.14) 


(dca) J(a;) 


(•••)g = 


(dw)JG(a;) . 
[{du)JG{uj) 


(4.17) 


In Ref. [52] it is shown that the solution of Eq. (I4.16p in the matrix form is, 

D = M(I + VM)“\ (4.18) 
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where I is the identity matrix, 

\\A^kfiq)dk,+k2,oh -q, 

2?tfci,fc2(q, -q, -Jy) = i^(Tpfci(q, z^)pfc2(-q, -^)) 

is the two-particle correlator, where the averaging is performed with the Hamiltonian 
of the system, 

(...) = lsp(e-^(^-^^)...)- (4.19) 

With this approximation the functional representation for Hint is, 

f^int =-^InSint =-^In JJ JJ JJ (^^z/fc(q)) ^ f(du) Joiuj)■ (4.20) 

S S q^O u k 

This Gaussian functional integral can be easily calculated. The result is, 

^int = -^In JJ JJ JJ (A^fc(q)) ' (l^^(q’ ' 


q^O k 


_ Lvvi 

2/3^^ "" ni^zc(q) 


(4.21) 


q^O 


From this expression it follows that it is necessary to know the effective inter-electron 
interaction for further calculation of Hint. Its Fourier-transform satishes the matrix 
equation (14.151) . This equation can be written as, 

^fci.fc 2 (q, = ^Z^fci(q)4i+fc2,0 + ^fci(q)£'-fci.fc(q, Z^)^fc,fc 2 (q, z^)- 

k 

Because, 


(4.22) 


L I L 

'2 2 


£/fci,fc 2 (q, z/, zi, Z 2 ), 


L _L 

■ 2 2 


(4.23) 


g{ci,u,zi,Z2) = X] e 

ki,k2 


— ifc222 


^fci,fc2(q, z^) 


and. 


L I L 

’2 2 


^fci.fc 2 (q> z/) = — dzi dz 2 e ‘*^"^2^(q, z/, Zi, Z 2 ), 


L —k 

' 2 2 


H(q, z/, 2 ;i,Z 2 ) = gifcizi-|-ifc2 22£)^ ^ 


k\,k2 


(4.24) 
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the equation 04.221) in (q, z)-representation has the form, 

5'(q, - Z 2 ) 


L I L 

'2 2 


SL^ 


dz / dz' z/(q, zi - z')D{ci, z/, z', z)g{(\, u, z, Z 2 ), 


(4.25) 


This integral equation was solved with different approximations in Refs. 

Consequently, the following expression for the thermodynamic potential 12 04.2p can 
be obtained, 

1 




(4.26) 


where. 


f!„ = -1 hiHo = - 15 ] 111 [1 + 


P 


p,Q: 


is the thermodynamical potential of non-interacting systeirU, 

1 det g{q,iy) 

= TTi 


q^O ly 


Ylr’^kiq) 


(4.27) 


(4.28) 


Therefore, the calculation of the thermodynamic potential in this approach requires 
knowledge of the effective inter-electron interaction with taking into account the presence 
of the dividing plane. 


4-2. The calculation of Qq 

Let US now calculate the thermodynamic potential of the non-interacting system, 

1 


12o = - - ^ In [1 -k 


(4.29) 


p,Q; 


Because here pi is the chemical potential of interacting electrons, this expression takes 
into account the Coulomb interaction indirectly. 

In order to perform the summation by p and a, we use the density of states 
calculated in the Appendix Appendix A Then, the thermodynamic potential is. 


Qo = -^ j dEp{E)ln [l + = 


1 SL \/2m3/2 


1 'j2rfV'l‘^d 


(3 2 

CO 

jdEPE\n[l + 


dE Ve In [1 + 

) 

dE In [1 + . 


m 




f Because the form of the thermodynamic potential flo as a function of g coincides with the 
thermodynamic potential of an ideal electron gas, the thermodynamic potential fio (14.2711 is called 
by us as ‘thermodynamic potential of non-interacting system’ though it indirectly takes into account 
the Coulomb interaction between electrons via the chemical potential g of interacting electrons. The 
same applies to the internal energy of non-interacting system Uq. 
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Integrating by parts each of the terms, we get, 


fin — 


OO 








2 V2m^/^d 


dE 


+ S- 


m 


dEE 


qP(e m) -|- 1 ' 47rh^ 

b 0 

In the limit of low temperatnres {fi —>■ cxo), we get the following expression. 


y(E-ii) _|_ ■ 


Hq f^0,bulk E 


O.surf; 


where. 


„ SLiV 2 m^E 

i‘0,bulk —-- 1 i-_9fc,q h — 


SL 


jPf 


(4.30) 


(4.31) 


2 ISvr^h^ 2 ISmvr^ 

is the extensive contribntion to the thermodynamic potential of the non-interacting 
system (it is proportional to the volume SL), which depends on the Fermi momentum 
Pp of interacting electrons and. 




0,surf 


= -5 


/ 4\/2m^/^(i 
\ ISvr^h^ 


-h 


5/2 _ 


m 

Snh? 




= S^^ 




TT 


dpF 


(4.32) 


is the surface contribution (it is proportional to the area of the dividing plane S). 

It should be noted that the thermodynamic potential of the non-interacting 
homogeneous system is 


o _ T/ 4\/2m^/2 5/2 _ 


ISvr^h^ 


ISmTT^ 


Pf, 


(4.33) 


which coincides with fl4.3ip . 


4-3. The calculation o/flint 


As shown in Eq. fl4.28p . the calculation of the thermodynamic potential requires 
to evaluate the determinant of the matrix effective inter-electron interaction. It is 
challenging, since the order of the matrix is inhnite. We can use the well-known identity 
(see, for example, [6T] ) 

In det A = Sp In A, 


with further expansions in series In A, next calculating the trace of each matrix term in 
the series and following summation of the series. 

To avoid this, we propose a different procedure (it is equivalent to the approach 
outlined above). The essence is to build a differential equation for the unknown quantity 
flint- To this end, we introduce the function instead of Jg{u) (14.141) . which 

depends on the parameter A: 


JG(a;, A) = exp 


fci,fc2 


(4.34) 
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that depends on the parameter A, 

= i-iL’^ki{q)y^hi+k2,o - AZlfci,fc2(q,z/), (4.35) 

moreover !)• 

Then, f2int and Sint will depend on this parameter as well, 

^int(-^) ^la.^ini,(A) 

= n n n A), (4.36) 

^ qT^O u k 

moreover 

^int ~ ^int(l)- 

We need to perform differentiation of flint (A) with respect to the parameter A, and 
arrive at the result, 

dflint(A) 1 d 1 1 dSint(A) 

dA “ /3dA /5Sint(A) dA 

=-;55-7^ n n n >< 

jidu)JG{uj,\) T)fcj,fc2(q,z^)a;fci(q,z^)a;fc2(-q,-z/) = 

cij^0,u 

ki,k2 

= DkiMi^^yykM,y^k2{-(i,-y)Gy)^ 

qf^0,iy 

k\,k2 

where we have introduced the average, 

f {du)JG{u},X)... 

{...)Gy) = ^r -• (4.37) 

J (dcn)JG(a;, A) 

Both averaging fl4.17l) and fl4.37|) coincide when in the latter the parameter A is 1, 
(•••)g = (...)g(1) 

Therefore, flint (A) satishes the differential equation of the hrst order, 

dfl^n^(A) = Y Dk^,k2i^^yyk^{q,yujkfi-q,-iy))G{y- (4.38) 

k\,k2 

In order to obtain an unambiguous solution of the differential equation of the hrst order 
it must be supplemented by the single additional condition. It is worth noting that, 

fiint(o)= -4 In n n n 

q^O k 
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X {du)exp -- W(q,z^)^^-fc(-q, -i^) 


C[:^O^U,k 


= -l„l = 0. 


(4.39) 


We easily find the solution of the Cauchy problem (I4.38p . (14.391) : 

1 

f^int = f^int(l) = ^ Zlfci,fc 2 (qw) [{(^kfiq,T^)uJkfi-(l,-iy))^{X)dX. (4.40) 


ci^O,^ 

ki,k 2 


Average of a;fc^(q, i/)a;fc 2 (—q, —u) with the Gaussian distribution X) yields, 

(a;fci(q, z/)a;fc2(-q,- i/))g(A) = ^^fci,fc2(qW, A). 

Thus, we have obtained a convenient expression for the calculation of Gjnt: 


Gint — 


i 

^fcnfc2(qw) y'^fci,fc2(qw, A)dA. 


(4.41) 


(4.42) 


q^0,v 

ki,k 2 


Using the relations fl4.23p (which holds for the effective inter-electron interaction 
dependent on A) and fl4.24p . we get. 




2 SL^ 


/ ^^2 iA(q, I/, zi, 2 : 2 ) g{q^,i',Zi,Z 2 ,X)dX. (4.43) 


q7^o,i/ ^ 


_ _L 
2 2 


Further evaluation of Gjnt using this formula should be carried out numerically. In 
order to obtain the analytical results, we make the following approximation: 

• D ^ 971°, namely we apply the random phase approximation; 

• g{q,i',Zi,Z 2 ,X) ~ 5 f(q, 0, Zi, ^ 2 , A) = g{q, Zi, Z 2 , X), i.e. we neglect the dependence 
of the effective inter-electron interaction on Bose frequency u. 

Then, the expression for Gjnt is simplihed, 

1 


L .L 

’2 2 




int 


25L2 


/du/d^2 5^9It° {ci,u,zi,Z 2 ) g{q,Zi,Z 2 ,X)dX. 
q^o L L n 


(4.44) 


In this expression the summation by the frequency z/ only applies to the function 


m^{q,u,Zi,Z2) = ^ 


^«i(p) - ^«2(p - q) 

^ p,tT;'a2 (p - q)' 

and can be performed analytically, 

^97l°(q, iy,Zi,Z 2 ) = -L^'^nafip)\(pai{zi)\‘^6{zi - Z 2 ) 

U P,fA:i 


</5ai(u)V?a2(u)V^^^2)V?„,(^2)(4.45) 


+L2 ^ n^Av)no,fiv-q)vlfzi)Lp^lzi)^llz2)ip^fz2). 


(4.46) 


P,Q1,Q2 
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The effective inter-electron interaction g{q^, zi, Z 2 , X) is the solution of the integral 
equation: 


g{c^, zi, Z2, X) = iyiq,zi - Z 2 ) 

+L/2 -\-Lj2 

/S 


(4.47) 


+ X dz dz' iy{q\zi - ^')a)l°(q, 0, z', z)g{ci, z, Z 2 , A). 

-L/2 -L/2 

Substituting Eq. fl4.46p into Eq. f|4.44p . and Eq. fl4.44p in Eq. fl4.26p . we hnd the 
thermodynamic potential, 

+f 1 

^ = ^0 - ^(A^)oX]i^(q,0) + jdz\(pa,{z)\^ jg{q,z,z,X)dX (4.48) 


q^O 


p,a 



2 2 


Taking into account the expressions for one- and two-particle distribution functions 
of electrons in the semi-inhnite jellium [50] . 

Fi{z) = iT’aWI V(p), (4.49) 


F^{r\\,z^,Z2)=F^{z^)F^{z2) 


S^N)l 


e'P’'iin„,(p')n„2(p' 


p,«i 

p',02 


(4-50) 


the thermodynamic potential can be represented as, 

+# 

2 

~ 2 ~ 2 - 
+ Jdr\\JdziJdz 2 {F^{r\\, zi, Z 2 ) - F^{zi)F^{z 2 )) J dX g{r\\, zi, Z 2 , X), (4.51) 

S _L 0 

2 2 

where, 

g{r\\,Zi,Z2,X) = i^e'‘^^ii^(q,2;i,2;2,A) (4.52) 

q 

is the effective inter-electron interaction in the coordinate representation, which depends 
on the parameter A. 

It should be noted that the expressions fl4.49l) and fl4.50p coincide by form with the 
expressions for the distribution functions of electrons without the Coulomb interaction, 
but these distribution functions depend on the chemical potential /i of interacting 
electrons. 


j dA^(r||,^,z,A)|r||=o 
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Thermodynamic potential in the case of the infinite barrier model. The form of the 
surface potential Kurf(^) must be specihed to perform further calculation of the 
thermodynamic potential according to the expressions (14.481) or (I4.5ip . We use the 
inhnite barrier model for the surface potential, namely, 


14urf(2;) = 


oo, z > d, 

0, z < d. 

The wave functions and the corresponding energy levels for the model are. 


iPa{z) = 


sin {a{d — z)), z ^ d, 
VL + 2d 1 0, z > d, 


S/y 


hfa‘^ 

2 m 


where. 


Tin 


a = 


-, n = 1, 2, 3,... 


(f + 

The one-particle distribution function of the model is, 

3cos {2pp{d — z)) 3 sin (2pF(d — ; 2 )) 


Ff{z) = 


1 + 


(4.53) 


(4.54) 


(4.55) 


9{d-z). (4.56) 


{2pF{d-z)f {2pF{d-z)f 

This expression for d = 0 coincides by form with the one-particle distribution function 
without Coulomb interaction [621 - I64] . but the expression fl4.56p takes into account the 
Coulomb interaction through the Fermi momentum pp of interacting electrons. 

Using the technique to solve the integral equation (I4.47p (see [53l|5l]), we obtain 
the following expression for the effective inter-electron interaction g{q^, zi, Z 2 , A), 


g{q,zi ^d,Z 2 ^ d, A) = 
g{q, Zi ^ d,Z 2 ^ d, A) = 
g{q,zi ^ d,Z 2 ^ d, A) = 
g{q,zi ^d,Z 2 ^ d, A) = 


27re^ 

^)L 

27re^ 


g-Q(A)pi-Z2| 


q 


Q{X) + q 
dvre^ 
Q{X) + q 


QW_q Q{X)(zi+Z2-2d) 

Q{X) + q 

-q\zi-Z2\ _ QjX) q q(zj^^z2-2d) 

Q(X) + q 

_ pQW(.Z 2 -d)-q{zi-d) 


fi{\){zi-d)-q(z 2 -d) 


where. 


Q{X) — + X XtfL(^), 

1 l-\-x 

Tix) = - H-^-In - 

2 Ax 1 — X 


^TF = Y inverse Thomas-Fermi radius of screening, and ob is the Bohr radius. 

After the summation by the momenta p in fId.dSp . we get, 

^ = ^bulk + ^surf, (4-57) 

where the hrst term is the extensive contribution to the thermodynamic potential (it is 
proportional to the volume of the system SL), the second term is the surface contribution 
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(it is proportional to the area of the dividing plane S). The extensive contribution to 
the thermodynamic potential is, 

f^bulk f^Ojbulk T '^f^bulk; (4.58) 

where, 

f^0,bulk ^5 f A |-n^ 

is the extensive contribution to the thermodynamic potential of the non-interacting 
system per unit volume (see the expression fl4.3ip h This contribution depends on the 
Fermi momentum pp of interacting electrons. Affbuik has the form, 

CO 1 


where. 


Afibuik _ ey-f 
SL/2 ~ 67r2 


dq 


dA 


Q 


Q(A) 


0 0 

CO CO CO 


27r4 


‘dqq /dai/da 2 J(g, «!, 02 ) d\ gfiq, ai, a 2 , X), (4.60) 


0 0 0 


TTCf, C2 > Cl, 
-,2 


0 ^q < |ci - C 2 I, 


J{q,ai,a2) = < 


TTC^, Cl > C2, 

/(ci, C 2 , q) + /(c 2 . Cl, q), |ci - C 2 I ^ g < Ci C 2 , 


0 , 


q> Ci-\- C 2 , 


Ci = \Pl-afi C2 = \lpl-al, 


c c‘^ 1 ^ arcoin ^'"^^^^' ^ - cj + / 2 (cf - cj-h q^fi 

f(ci,C 2 ,q) Cil^ arcsm J J c, 


“i) “2, A) 


Q^(A) -f -f 

(Q^(A) + aj + al)'^ - Aafaj 


It should be noted that the dividing plane has no effect on the expression given 
by Eq. (14.581) and actually this expression is the thermodynamic potential of the 
homogeneous system per unit volume. However, the thermodynamic potential of the 
homogeneous system can be calculated in a similar manner, based on the Hamiltonian 
of the inhnite jellium fl2.5p . These calculations are much simpler and at the similar level 
of approximations we obtain, 

_ 00 1 

j dA5(unif(g, A) - i2{q) 

0 0 
00 1 

” 12 ^ /d??^'^nif(g) ydA gunif(g, A), (4.61) 

0 0 


a 


unif 




V 


IdmTT^ 


Pf + 


Pf 


dq (fi 
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where, 


^unif(g) = 




0 , 


q ^ 2 pF, 


^{q) = is the three-dimensional Fourier-transform of the Coulomb interaction, 




47re^ 


-|- A x^pL(2pp) 

is the three-dimensional Fourier-transform of the effective inter-electron interaction of 
the homogeneous system, that depends on the parameter A. 

Integration by the parameter A in the expression given by Eq. (14.611) can be easily 
performed and as a result we obtain. 


n 


unif 






V ISmvr^^^ Svr^ 


dg 


^ _ 9^(1 - ^nif(9)) 




q 


2pF ' 


(4.62) 


Although the expressions fl4.58p - fl4.60p and fl4.62p are different by form, they both lead 
to the same result, namely to the thermodynamic potential of the homogeneous system 
in the random phase approximation. 

The surface contribution to the thermodynamic potential has the form. 


Alsurf All 


where. 


n 


0,surf 


0,surf 


hYF 


miT^ 


AQ. 


surf ? 


71 dpF 


(4.63) 


(4.64) 


is the surface contribution to the thermodynamic potential of non-interacting system 
per unit area (see expression fl4.32l) ). which depends on the Fermi momentum pf of 
interacting electrons, 

CO 1 


An 


surf 




127r2 


e^dpl 


dg 


dA 


q 


Q(A) 


dg /dA 


0 0 

g Q(A)-g 


+ 


e^PF 

167r 


dg 


0 

CO 




Q2(A)g(A) + g L 

CO 1 

dqqj daijda 2 J{q,ai,a 2 ) Jd\ g 2 {q, ai, a 2 , X), 
0 0 0 

g(A)-g 


l + ^-^^iM^i±^arctan 


^Pi 


QW 




92 ( 9 , ai, a 2 , A) = |2gi(g, ai, 03 , A) d + ^ ^ 

The parameter d for the inhnite barrier model is, 
dvr 


(4.65) 


16g(A)a^a^ 


[(g2(A) + a\ + a. 


2\2 


Aafalf 


d = 


8 pf' 


(4.66) 
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It is worth mentioning that the second term in the expression given by Eq. fl4.65p 
(it is linear on the chemical potential and is quadratic on the Fermi momentum p-p) 
contains the divergent integral. However, as we will see below, this divergent integral 
disappears in the calculation of the internal energy. 


5. Internal energy 
5.1. General expressions 

By using thermodynamic potential and the Gibbs-Helmholtz equation generalized for 
the case of variable number of particles, we obtain the internal energy of the system U, 


^ ^ ^89 ^dpi' 


(5.1) 


At low temperatures 0 —)■ 0, the second term of the r.h.s. of this equation vanishes and 
we get, 

U = Vt + pi{N), (5.2) 


where we have used the relation 


{N) = -Sp (e^(^-^^)iV) = 

S ^ ^ dpt. 


(5.3) 


According to the Eq. (I4.57p . thermodynamic potential can be divided into the 
extensive and surface contributions. Then we get. 


where. 


/\T\ _ ^(^bulk + hlsurf) _ jyj- jyj- 

\tv; —- — -— iVbuik + JVsurf, 

^f^bulk 


^bulk — 
-^surf 


dp. 

cAlsurf 


and 


dp ’ 

U = Hbulk + f4urf, 

fdjulk f^bulk P f^bulk T P-hlhulk: 

dp 

Usuvi ~ f^surf h o ~ f^surf T p-hlsurl 
dp 


(5.4) 

(5.5) 

(5.6) 


(5.7) 

(5.8) 


are the extensive and surface contributions to the internal energy, respectively. 

The chemical potential p is the solution of the equation (15.311 . By using (15.411 in 
the thermodynamic limit, we get 


lim 


(TV) 


= lim 


Ny 


bulk 


3 1 


n,s,l^ooSL/2 n,s,l^oo SL/2 dvrrf ’ 
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where r^ is the Wigner-Seitz radius in units of the Bohr radius ob- Moreover, in the 
thermodynamic limit the summand does not affect the chemical potential fi (but 
affects the surface contribution to the internal energy t4urf) and the equation for pi can 
be presented as, 

W = (5.9) 

According to Eqs. fl5.7p and flS.Sp . in order to calculate the extensive t/buik and the 
surface f4urf contributions to the internal energy we need to evaluate the extensive A^buik 
and surface Nguri contributions to the average of the number operator of electrons (iV). 


5.2. Average of the number operator of electrons and the chemical potential 


According to Eqs. (15.51) and (15.6|) . A'buik and Asurf can be calculated by taking the 
derivatives of hlbuik and hlsurf with respect to the chemical potential /r, respectively. 
However, it is easier to use the functional representation of the thermodynamic 
potential hi (see (I4.20p i and compute the derivative of this expression (I4.26P with respect 
to the chemical potential pi. As a result, we hnd that. 


1 


(Af) = Wo + ^ 5^ Kq. 0) 


q^O 


d{N)o ^ 1 aSint 

dfJ. /3^int 


9/x 




q^O 


H' 'int — 1^ 




= {N) 


q^^O u k 

1 


dp, 


UJkfici, u)ukfi-q,-n) 


k\,k 2 


0 


2^ 






q^O 


dp 2(3 ^ dp 

^^ 0 ,u 

ki,k 2 


{Ukfiq, u)ukfi-q, (5.10) 


where averaging (.. .)g is performed according to Eq. (14.171) . 
Considering that, 

(a^fci(q, iy)ujkfi-ei, -i^))g = 


the expression (I5.10p can be rewritten as. 


q^O 


dp 


2S 


dp 


9ii,fa(q. ")■ 


(5.11) 


(5.12) 


ki,k2 


Taking into account Eqs. (I4.23P and (I4.24p . we obtain. 


1 
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25L2 


~ 2 2 




_ _L 
2 2 


gZl(q, z/, ^ 1 ,^ 2 ) 
dp. 


g{q,u,Zi,Z2), (5.13) 


In order to simplify this expression we make similar approximations, as in the calculation 
of the thermodynamic potential, namely: 

• 5'(q, Zi, Z2) ^ 5'(q, 0 , Zi, Z2) = fi'(q, Zi, Z2). 

Then the expression given by Eq. fl5.13|) is greatly simplihed and. 


qT^o 


dp 


-L ,L 
2 ~ 2 


25L2 


/dzi /dz2 

qT^o \ 


k —k 

2 2 


gajl°(q, z/, ^ 1 , 2 : 2 ) 
dp 


g(q,Zi,Z2), (5.14) 


where the summation by the frequency n applies only to the derivative of the two-particle 
correlator with respect to the chemical potential p and can be performed analytically. 




991l°(q, u, Zi,Z2) 


dp 




PjOil 


dp 




P,Ol,Q2 


dng^ (p) 
dp 


(n„ 2 (p - q) + Ugfip + q))^lizi)ip^fzi)^plfz 2 )iPgXz 2 ). (5.15) 


Substituting Eq. fl5.15|) into Eq. fl5.14p . we hnd, 

W = (iV>0 + T „(q, 1 ^ ^ S"»(P) 


q^O 


dp 2S dp 

^ qT^O p,a ^ 


dz\ipg{z)\^g{ci,z,z) (5.16) 


1 ^ ^ d^fip)ngfip - q)) 


L I L 

’2 "*”2 


2^ 


q^O p,01,02 


dp 


dzi / d ^2 ^Uzi)^Jzi)^*Jz 2 )^gfz 2 )g{q, zi, Z 2 ). 


k —k 

■ 2 2 


Taking into account the expressions for the one- fl4.49p and two-particle fl4.50l) 
distribution functions of electrons in the semi-infinite jellium [50], we get. 


w = Wo + ^ E "(i'°) 


d{N)„ IS f, d((JV)oF«(z)) 


q^O 


dp 2 V 


dz 


dp 


■g{r\\,z,z)\ ^ 


1 5 

'2E2 


. d{{N)liF^{v\\,z,,Z2)-F,%z^)Ff{z2))) , 

drii / dzi I dz2--5'(r||, 2:1, Z 2 ), 


s 

2 2 


dp 


(6.17) 


where g{Y\\,zi,Z 2 ) is the effective inter-electron interaction in the coordinate 
representation (note that the transition to the coordinate representation is similar to 
Eq. (1132D). 
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5.3. Calculation of {N)q andd{N)Q/dpL 


By using the thermodynamic potential of the non-interacting system hlo (see Eq. (I4.3n|l ). 
we calculate {N)o at low temperatures, 




^(f^0,bulk “1“ f^0,surf) 


where, 




dfi dp, 

^hlobuik SL2y/2m?C 


= N: 


O.bulk 


+ N^ 


0,surf 5 


0,bulk 


dp 2 

is the extensive contribution {N)q, and 


-p3/2 = 


SL 

2 37r2 


N, 


dfl, 


0,surf — 


0,surf 


dp 


= S 


2 \/ 2 m^^‘^d 


-p 


3/2 


m 


Anh? 


P\=S^-^ 


Pp f dpp 71 


TT^ 


(5.18) 


(5.19) 


is the surface contribution {N)q. The latter result reduces to the one obtained previously 
in Refs. [HSIIZI], if we put d = 0 in Eq. fl5.19p . 

By using {N)q, we calculate d{N)Q/dp 

d{N)o dNp 

,bulk dNp 

,surf 

dp dp dp 

where the extensive contribution is, 

^-^o.buik SL \f2rrdl‘^ SL m p-p 

dp 2 ^ 2 ^2 TT^ ’ 

and the surface contribution is. 


(5.20) 


dN, 


0,surf 


dp 


= S 


[ y/ 2 wfi/‘^d 
\ 7r2h,3 




1/2 _ Jff_ 

Anh^ I 


= S 


m 


(^dpp - 


(5.21) 


It should be noted, that the following relations are valid for the non-interacting 
homogeneous system. 


^0,unif — y 


pI 

37r2’ 


(5.22) 


^A^0,unif ^ yfifiPF 

dp h2 7r2 ’ 

They coincide with the extensive contributions given by Eqs. 
respectively. 


(5.23) 

flHTSll and flT^ . 


Average of the number operator of electrons for the infinite harrier model. For further 
calculation of fl5.16l) or fl5.17p . as in above, we consider the inhnite barrier model fl4.53l) 
of the surface potential 14urf(2^)- After the summation by momenta p in Eq. (I5.16p . we 
get, 

Abuik = Atg^buik + AA^buik, 


( 5 . 24 ) 
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where, 


ATf 


0,bulk 


pI 


SL/2 3.^ 

is the extensive contribution to the average of the number operator (N) of the non¬ 
interacting system per unit volume (see expression fIS.lSp b This contribution depends 
on the Fermi momentum pp of the interacting electrons. AiVbuik has the form, 

OD 

AAbuik Pf f, L q 

SLI2 “27r2aBy n Q 

0 


TT^Ob 


dqq dai / das I-{q, «!, 02 ) gi{q, «!, ^ 2 ,1), (5.26) 




0 , 


arccos 


/_(g,ai,a2) = < 


0 , 




arccos 


0 , 


g^+cf-c| 
2cig ’ 


g^+cf-c| 
2cig ’ 


0 ^ ^ Cl - C2, 

Cl - C2 < g ^ Cl -h C2,; 
Cl + C2 < g < 00, 

0 ^ g ^ C 2 - Cl, 

C 2 — Cl < g ^ Cl -|- C 2 ,J 
Cl + C 2 < g < CO, 


Cl > C2, 


C2 ^ Cl, 


(6.27) 


where. 


Cl = 


pI 


at 


'1) C2 — Y FF ~ '^2- 

It should be noted that the dividing plane has no effect on the expressions fl5.24p - 
fl5.26p and therefore Eq. fl5.24p is the number of electrons of the homogeneous system per 
unit volume. However, the number of electrons of the homogeneous system can similarly 
be calculated from the Hamiltonian of the inhnite jellium (12.5p . These calculations are 
much simpler and with similar approximations we obtain. 


pI 


Q^n 


N, 


unif 


Pf 


+ 


Pf m 


where. 


dgg^ (iy(q) -gunif(g)) + 


Pf m 


dq (5.28) 


'unif 


(q) = 


^unif(g) = 


2 , 

1 - 

2pp’ 

0 , 

dvre^ 


^-1, 


2pF 

-1 < ^ ^ 1 , 

2pF ’ 

- 2 - > 1 
2 PF ’ 


g2 -1- x|pL(2pp) 

is the three-dimensional Fourier-transform of the effective inter-electron interaction of 
the homogeneous system. 
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Substituting the expressions for the three-dimensional Fourier-transforms of the 
Coulomb interaction i/(g) and the effective inter-electron interaction gunnid) into Eq. 
(.I5.28P . we get, 

Nnnil _ Pf PF 7 (^TfL(^) + 

The expressions fl5.24p - fl5.26p and fl5.29l) are nonlinear algebraic equations for the 
chemical potential /i (/i = Although their form is different but they yield the 

same result, because the dividing plane does not affect the chemical potential. 



Figure 1. The chemical potential of the interacting electrons y (solid line) and the 
chemical potential of the non-interacting electrons /tq (dashed line) as a function of 
the Wigner-Seitz radius. 


In Figure[I]the chemical potential as the function of the Wigner-Seitz radius rg is 
presented. The nomenclature of lines is given in the hgure caption. It is worth noting 
that the chemical potential is the solution of nonlinear algebraic equations fl5.29p . and 
the chemical potential of non-interacting electrons is, 

'97rV/^ 1 

— ' -%• 


ho — 


4 7 (r-s/aB)^ 


(5.30) 


It can be seen that taking into account the Coulomb interaction leads to a decrease 
of the chemical potential of electrons. 

In Figure|2]the parameter, d, d = 37r/(8pF) as a function of Wigner-Seitz radius 
Tg is given. The parameter d is the distance from the inhnite potential barrier to the 
dividing plane {z = 0). The conclusion from the curves is that taking into account 
the Coulomb interaction between electrons leads to an increase of this distance and its 
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Figure 2. The parameter d as a function of Wigner-Seitz radius (solid line is for 
interacting electrons whereas the dashed line is for non-interacting electrons). 


nonlinear dependence on rg, whereas the parameter d for the non-interacting system is 
a linear function of rg. 



Figure 3. The one-particle distribution function of electrons as a function of the 
electron coordinate normal to the dividing plane at rg = 2 cb (the solid line is for 
interacting electrons whereas the dashed line is for non-interacting electrons). 

In Figures[3] and HI the one-particle distribution function of electrons fl4.56p is 
presented for the following values of Wigner-Seitz radius: Ug = 2aB and Ug = hae- The 
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Figure 4. The one-particle distribution function of electrons as a function of the 
electron coordinate normal to the dividing plane at rg = 6 ob (the solid line is for 
interacting electrons whereas the dashed line is for non-interacting electrons). 


one-particle function of electrons (solid lines) depends on the chemical potential which 
is the solution of nonlinear algebraic equations fl5.29p (here the Coulomb interaction is 
taken into account). The positive charge is located at z ^ 0. It can be concluded that 
taking into account the Coulomb interaction leads to an increase of period of damped 
oscillations of the one-particle distribution function around its value in the body of the 
metal which equals to unity. 

The surface contribution to the average number of electrons (N) has the form, 


^surf — 
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( 5 . 31 ) 


( 6 . 32 ) 


is the surface contribution to the average number of the non-interacting electrons 
system per unit area (see Eq. (I5.19p ). that depends on the Fermi momentum pp of 
the interacting electrons, and 
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It should be noted that the second term in Eq. fl5.33l) . (which is linear on the 
chemical potential and is quadratic on the Fermi momentum pp) contains the divergent 
integral. However, as we will see below, this divergent integral disappears in the 
calculation of the internal energy. 

Also it should be noted that if we substitute the value of d (I4.66p into the expression 
fl5.32p . it vanishes. If we take into account the Coulomb interaction between electrons, 
then iVo,surf 7^ 0. 


5 . 4 . Internal energy and surface energy 


Substituting Eqs. fl4.5ip and (15.17p into Eq. (15.2p . we obtain the internal energy. This 
expression can be considered as one possible energy functional that in contrast to the 
functionals used in the density functional theory, depends not only on the one-particle 
distribution function of the electrons, but also on the two-particle distribution function 
and the effective inter-electron interaction. 

The extensive contribution to the internal energy per unit volume follows from the 
substitution of Hbuik and A^buik (or Hunif and A^unif) into Eq. (15.7p . Then we get, 
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where the hrst term of the expression (15.341) (or (I5.35p ) is the internal energy of non¬ 
interacting system. 

By using the expressions for flsurf and Nsuri in Eq. (15. 8 p we obtain the surface 
contribution to the internal energy per unit area. We are interested in the case of low 
temperatures (6* —)■ 0). Then according to [67l[68], the ratio U^nri/S is the free surface 
energy a, and the magnitude of f/surf = <^3 is the work that is necessary for irreversible 
process of creating a new free surface S. The quantity a describes excess energy of 
surface area compared with the energy inside the body of the metal. Then, the surface 
energy a can be presented as. 


a = do -h Act, 


(5.36) 
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where, 


ao = 
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(6.37) 
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(5.38) 


By placing the parameter d fl4.66p in Eq. fl5.37p . we obtain, 


(To = 


hfpp 


(5.39) 


IGOmvr 

It coincides by form with the surface energy of non-interacting system 
but in Eq. fl5.39p . pp is the Fermi momentum that takes into account the Coulomb 
interaction between electrons. 

In Figured the dependence of the surface energy on Wigner-Seitz radius rg is 
presented. The solid line is the surface energy calculated using the formulas fl5.36p . 
fl5.38p . fl5.39p the solution of nonlinear algebraic equation fl5.29p yields the chemical 
potential. The dashed line is the surface energy of non-interacting system fl5.39p . the 
dash-dotted line is the result of Lang and Kohn [5] . 

The result given in this hgure show that taking into account the Coulomb 
interaction leads to an increase of the surface energy compared to the surface energy of 
non-interacting system. In addition, the surface energy calculated by us is positive in 
the entire region of Vg and in the interval rg = 5.5 — Gap it coincides with the surface 
energy calculated by Lang and Kohn. 


6. Conclusions 

The general expression for the thermodynamic potential of the semi-inhnite jellium 
model is obtained by using the method of functional integration. The knowledge of the 
two-particle correlation function of electrons and the effective inter-electron is required 
for the practical application of this expression. 
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Figure 5. The surface energy as a function of Wigner-Seitz radius (the solid line 
is for interacting system, the dashed line is for noninteracting system whereas the 
dash-dotted line is the result of Lang and Kohn 0)- 


It is shown that taking into acconnt the Coulomb interaction between electrons leads 
to a decrease of the chemical potential. It is also shown that the surface contribution 
to the thermodynamic potential does not affect the chemical potential. 

By using the inhnite barrier model, the extensive and surface contributions to 
the thermodynamic potential, the average of the number operator of electrons and the 
internal energy are obtained and studied at low temperatures. 

The influence of the Coulomb interactions between electrons on the behavior of 
the one-particle distribution function is studied as well. We obtained that taking into 
account the Coulomb interaction between electrons leads to an increase of the period of 
damped oscillations around its average value in the body of the metal. 

It is found that taking into account the Coulomb interaction between electrons, 
the distance between the dividing plane and the surface potential as a function of the 
Wigner-Seitz radius looses its linear behavior. Namely, it grows faster. 

Based on the expression for the surface contribution to the internal energy and 
modeling the surface potential by the inhnite barrier, the surface energy is calculated 
at low temperatures. It is shown that taking into account the Coulomb interaction 
between electrons results in a growth of the surface energy. The surface energy is 
positivein the entire range of Wigner-Seitz radius. In the interval Vg > Sob, the surface 
energy calculated by us is in a good agreement with calculations of Lang and Kohn [5] . 

Our calculations of the chemical potential and surface energy demonstrate that 
taking into account the Coulomb interaction between electrons is very important in 
the region of small rg, and that the inhuence of Coulomb interaction on these values 
decreases with increasing rg. 
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Appendix A. Density of states 


Let us calculate the density of states of electrons without Coulomb interaction, 

p(i!)= (A.l) 


P,q: 


where according to fl3.2p 


Eaip) = 


2m 


In the thermodynamic limit (S' —)■ cxo and L —)■ oo, the sum can be replaced by the 
integral, according to the Euler-Maclaurin formula 
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where Bk are the Bernoulli numbers. 

We perform summation over two-dimensional vector according to Eq. (13.dh and 
obtain. 
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(A.2) 


where two possible orientations of the electron spin is taken into account. 
The summation over three-dimensional vector p yields, 
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Now we consider the summation over a. Eq. f|4.55p implies that. 
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The transition from the sum over p to the integral is performed according to Eq. 
flA.2p . Then the density of states flA.ip is, 
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Transformation from the snm over a to the integral according to Eq. flA.4p leads to, 



0 



(A.5) 


where the hrst term (which is proportional to the volnme SL) is the extensive 
contribntion, and the second term is the snrface contribntion to the density of states. 
This expression coincides with the expression for the density of states, obtained in 
Ref. [2B], if we pnt d = 0 in Eq. flA.Sp . 

It shonld be noted that the density of states of non-interacting homogeneons system 


IS, 


Punif(E) = 


which coincides with the hrst term of Eq. flA.Sp . 
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